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We study the effects of diffusion on a A-gradient echo memory, which is a coherent optical quantum 
memory using thermal gases. The efficiency of this memory is high for short storage time, but 
decreases exponentially due to decoherence as the storage time is increased. We study the effects 
of both longitudinal and transverse diffusion in this memory system, and give both analytical and 
numerical results that are in good agreement. Our results show that diffusion has a significant effect 
on the efficiency. Further, we suggest ways to reduce these effects to improve storage efficiency. 



I. INTRODUCTION 

Quantum memory is an important tool in many quan- 
tum information protocols, including quantum repeaters 
for long-distance quantum communication pQ, and iden- 
tity quantum gates in quantum computation [2 . Numer- 
ous optical quantum memories have been developed, in- 
cluding electromagnetically-induced transparency (EIT) 
based quantum memory [3j 0] , far-detuned Raman pro- 
cess memory [5 , 6 , and photon-echo quantum memories: 
controlled reversible inhomogeneous broadening (CRIB) 
memory [3 [8], atomic frequency combs (AFC) memory 
[9], and gradient echo memory (GEM) [T0HT2] . A review 
of these schemes can be found in [13]. Of these schemes, 
the most impressive efficiency so far attained experimen- 
tally is 87% by A-GEM scheme [14] using warm rubid- 
ium vapor. In this paper, we will examine the effects of 
atomic diffusion on the A-GEM system, which may limit 
this efficiency for larger storage times. 

A-GEM is a memory using a 3-level, A-type atom (Fig. 
[I]). The input optical pulse couples the two metastable 
lower states through a control field. The excited state is 
coupled in the far detuned region, so the 3-level atoms 
can be treated as effective 2-level atoms. These effective 
two-level atoms have linearly increasing atomic Zeeman 
shifts along the length of the storage medium. The pulse 
is first absorbed, then by simply reversing the sign of the 
magnetic field, the pulse is retrieved in the forward direc- 
tion. The incident signal field is converted to a collective 
atomic excitation known as a spin wave, which is dis- 
tributed as a function of position. The Brownian motion 
of the gaseous atoms will cause diffusion, which will dis- 
turb spatial coherence of the atomic spin- wave, leading 
to decoherence. For a plane wave, only axial diffusion is 
important, but transverse diffusion becomes significant 
when a realistic beam profile is included. There has been 
recent interest in the effects of diffusion in the EIT [T5l - 
H7j . In this work, we study the effects of diffusion in 
A-GEM system, giving analytical and numerical results. 
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FIG. 1. Level structure of A-type 3-level atom. 

We also suggest ways to reduce these effects to improve 
storage efficiency. 

II. A GRADIENT ECHO MEMORY 

We consider a medium consisting of A-type 3-level 
atoms with two metastable lower states as shown in 
Fig. [I] The ground state |1) and the excited state |3) are 
coupled by a weak optical field, the positive frequency 
component of the electric field is described by the slowly 
varying operator 

E(r, t) = J2 ]J^a^t)e ih - r e- ikoZ e iuJot (1) 

k 

with detuning A, where V is the quantization volume, 
Uq is the carrier frequency of the quantum field and 
ko = uo/c. The excited state |3) is also coupled to the 
metastable state |2) via a coherent control field with Rabi 
frequency Q c and a two photon detuning S. This two pho- 
ton detuning is spatially varied 5(z,t) = r](t)z, with time 
dependent gradient rj(t). Then the interaction Hamil- 
tonian in the rotating frame with respect to the field 
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frequencies is 

H = Y^[H^33 + M(z„,t) ( x. ( ™ ) 
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where = P^J 2 nt v * s tne atom-field coupling constant 
with p being the dipole moment of the 1-3 transition, and 
crffi = \p) n {v\ is an operator acting on the n-th atom at 
r n = (xmUn, z n ). We assume that initially all atoms 
are in their ground state |1). We transform to collective 
operators, which are averages over atomic operators over 
a small volume centered at r containing N r ^> 1 particles, 
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From the Heisenberg-Langevin equations in the weak 
probe region (an ~ l,c"22 — ^"33 — 0), we get the 
Maxwell-Bloch equations [18], 
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where 7^ are the decay rates, g = p 



77^-, and N is the 

2he Q ' 



atomic density. We have assumed that g^ 



P 
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and fi c (r) = Q c e 1 cZ . We also omit the Langevin noise 
operators since here we are more interested in the deco- 
herence caused by diffusion. This is equivalent to making 
a semiclassical approximation for the electric field and 
the atomic_coherences. 



In Eq. (4) 
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is the diffraction term, and gener- 



ally, the diffraction effects can be neglected [19]. Notice 
that we are here considering the regime t p ^> L/c, where 
2t p is the temporal width of the signal, and 2L is the 
length of the medium. This allows us to neglect tempo- 
ral retardation effects, i.e., we can neglect the temporal 
derivative in the third equation of Eq. Q. Also, since the 
atoms are far detuned ( A ^> 713, O c ), we adiabatically 

eliminate the fast oscillations and set &^ = 0. Then we 
have <Ji3 = (ge zk ° z E + Vt c e lkcZ G\^) /A, and we get the 
reduced Maxwell-Bloch equations, 
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Here we neglect decay, i.e. 712 —> 0, since we consider the 
storage time much less than I/712. 



III. DIFFUSION 

We now consider the effects of diffusion on the atomic 
state. In order to isolate the motional effects of diffusion 
from collisional dephasing, we assume that the collisions 
between atoms do not change the state of the atom. Then 
we derive the diffusion equation for the atomic density 
matrix p. Space is divided into volume elements with 
length Ar and center r. We associate a density matrix 
p(r, t) with atoms in this volume element, given by 

N r 



N r . 

J = l 

where N r is the atom number in volume centered at r. 
The total density matrix for the entire system is assumed 
to be the tensor product of these local density matrices. 

Diffusion causes an exchange of atoms between adja- 
cent volumes. During a short time At, a fraction e of the 
atoms in slice r migrate into slice r ± Ar. There is also 
atomic flux back into slice r from r ± Ar. We assume 
that the total number density of the atoms is uniform, 
so the state at r and t + At is described by the new den- 
sity matrix which is the average of the density matrix of 
atoms remaining in the volume and those that have mi- 
grated in to it. The diffusive component of the evolution 
is therefore 

p(r,t + At)=(l-2e)p(r,t) 

+ e(p(r + Ar, t) + p(r - Ar, t)) (6) 
=>dtp(r,t) =DV 2 p(r,t) 

where D = eAr 2 /At is the diffusion coefficient. With 
the same consideration, we get the diffusive component 
evolution for the atomic correlation functions 



(7) 



Now we introduce the interaction with optical fields. 
Since diffusion is caused by Brownian motion, this will 
lead to Doppler shifts in the various detunings. We 
now consider the interaction between the optical field 
and a single atom, and quantify the effects of these 
Doppler shifts. The atom moves at some random ve- 
locity, and there will be a Doppler shift for both the 
signal and control fields. So the detunings in Eq. Q be- 
come A = A + A Dopp , and S = S + S Dopp , with A , S 
the detunings for stationary atoms and Ad opp , Spopp the 
Doppler shifts. Typically, the one photon Doppler shift 
^Dopp <C Ao, state |3) is still far detuned. So the adia- 
batic elimination is still valid in the presence of the Brow- 
nian motion induced Doppler shift, and we can still re- 
duce the 3-level atom to an effective 2-level atom. The 
Maxwell-Bloch equation will still reduce to Eq. ([5|, but 
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with one photon detuning A = A + A d opp and two pho- 
ton detuning 5 = So + 5dopp- So, for the reduced two level 
atomic system, the diffusive Maxwell- B loch equation for 
the collective correlation (712(2, t) averaged over atoms in 
each volume is 



iS(z,t)a 12 (r,t) + DV 2 a 12 (r,t), 



oz cA 



<Ji2 (r,t) 



(8) 



We have absorbed the Stark shift into the two-photon 
detuning. Here our diffusive Maxwell-Bloch equation is 
consistent with the result in the EIT system p~5j [17] . 

Notice that the signal and control fields are co- 
propagating, so the Doppler broadening width for S is 
typically 1kHz, which is much smaller than the frequency 
width of the signal field (~ 1MHz), so we neglect this 
two-photon Doppler broadening 5d opp and replace 5(z,t) 
by 5 (z,i) = rj(t)z,z G [— L, L). 

For the one photon detuning A = Aq + A]j opp , after 
we make the adiabatic elimination, it will appear in the 
denominator (see Eq. pi), so 
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A]j opp \ 



The term linear in A]j opp will vanish when we average 
over many atoms in a volume centred at r, so we can 
replace A by Aq in our Maxwell-Bloch equation, with 
second order accuracy [typically (Ad opp / Aq) 2 ~ 10 -3 ]. 



A. Longitudinal diffusion 

For a uniform plane wave, transverse diffusion is irrel- 
evant. We replace r by z in Eq. ([8| and consider the 
longitude diffusion in a 1-dimensional model. Now the 
Maxwell-Bloch equation is 



( k °- kc ^ z E(z,t) 
i5(z,t)a 12 (z,t) + D\7 2 z a 12 (z,t), 



^ N W( A 
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(10) 
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We now investigate the longitude diffusion effects dur- 
ing the write-in process, the hold time and the read-out 
processes separately. 



To compute f ou t(t), we evolve Eq. (10) using r\ as in 
Fig. [5] (a). Following the method given in [10 , we first 
propagate E(z,t) and a 12 {z,t) forward with boundary 
condition E(z = —L,t < 0) = fi n {t) to find their values 
at time tn- Then we propagate E and cr 12 backward to 
time tui with final condition E(z — L,t > tjj) — fout(t)i 
and solve for f ou t(t) by matching the two solutions at 
time tu- 

Write: Consider the diffusion effects during the write- 
in process, we find that (see Appendix A) 



fout(tH + 1) 



ew(*?-(*'-*> 3 >/ in (-t)G 



(11) 



where ki = ^^+ko—k c — j; is the initial spatial frequency 
of 012(2, t), and 



IV. ANALYTIC CALCULATION AND 
NUMERICAL SIMULATION 

To quantify the effects of diffusion, we solve for the 
atomic dynamics. There are three distinct phases during 
the storage: write-in — to < t < 0, during which the signal 
is absorbed by the memory; hold < t < tn, during 
which the information is stored in the memory and the 
gradient is turned off; read-out tjj < t < tjj + to, during 
which the signal is emitted by turning on the flipped 
gradient. 

We quantify the effects of diffusion by the read-out 
efficiency e defined to be 



f t H H +t °\fout{t)?dt 
1° \fin{t)?dt 



where f ou t(t) = E(z = L,t > tu) is the output field 
and fi n (t) = E(z = —L,t < 0) is the input field. We 
solve for f ou t(t) both numerically and analytically, and 
consider the effects of diffusion in axial (longitude) and 
radial (transverse) directions separately. 
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is a phase factor, with g e ff = gQ c / A and j3 — 

For a pulse with Gaussian temporal profile f 7 
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we find 



J —to 



(12) 



Typically, Drj t p is very small, then the efficiency is 
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and rw 
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are dimensionless pa- 



rameters. For typical experimental parameters, aw — 1, 
then 
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(14) 
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We also numerically solve Eq. (|10| with diffusion during 
the write-in process, using XMDS [20 . We calculate the 
efficiency for different values of the diffusion rate D, input 
time ti n etc. The results are shown in Fig. [2j (points are 
numerical results, and the curve is Eq. (|l4|)). We plot the 
efficiency sw with respect to the rescaled dimensionless 
parameter r^y, so all the points with different parameters 
collapse on a single curve. 

Hold: During the storage time [0, tn], we find (see Ap- 
pendix A) 

fout(t H +t) = e- Dt »^-^ 2 f in (-t)G (15) 

For the above Gaussian shape input, the efficiency is 
given by 

e H = V^- 2 ™, (16) 

where an — p-/ + DtnV 2 ^ an d th = Dtnk 2 H are di- 
mensionless parameters, with kn = h — vUn- For typical 
experimental parameters, an — 1 and we have 

-2th 
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We also numerically solve Eq. ( 10 ) with diffusion during 
hold time, using XMDS. We calculate the efficiency for 
different values of the diffusion rate D, storage time tn 
etc. The results are shown in Fig. [3] (points are numerical 
results, and curve is Eq. (17)). We plot the efficiency sh 
with respect to the rescaled dimensionless parameter r#, 
so all the points with different parameters collapse on a 
single curve. 

Read: the diffusion effects during the read-out process 
are the same as the diffusion effects of the write-in process 
(see the appendix A), so we simply have 
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FIG. 2. The efficiency decay with respect to the dimensionless 
parameter tw for longitude diffusion during the w rite -in time, the 
points are numerical results, and the curve is Eq. ( [14} . 




B. Transverse diffusion 



FIG. 3. The efficiency decay with respect to the dimensionless 
parameter th for longitude diffusion durin g ho ld time, the points 
are numerical results, and the curve is Eq. jT7|). 



We now quantify the effects of diffusion for a beam 
with realistic transverse Gaussian profile. The efficiency 
for a 3-dimensional model is defined as 



J \fout(x,y,t H +t)\ 2 dxdydt 
J \f in (x,y,t)\ 2 dxdydt 



(18) 



Eq. (|8| can be solved in Fourier space k x ,k y , and also 
notice that 



J \fout(k x ,ky,t H +t)\ 2 dk x dk y dt 
f \fin(k , x-> i)\ 2 dk x dkydt 



(19) 



Eq. ([8]) can be reduced to a quasi- ID problem, and 
can be solved as before (see Appendix B). For transverse 
diffusion, the output pulse will be 



t H + t) 



6 Tfe g 7fc fin(yk x ^ky^ t)Cr. 



(20) 



where j k = D(k% + k 2 ). 

If the input pulse has both Gaussian temporal and 
transverse profile, 

f m (x,y,t) = Ae-^+y^e-^ 2 ^, 

then jktp ~ Dtp /a 2 , which is typically small. Thus the 
memory efficiency is 



1 + TL 



(21) 



where r±_ = 4:D(tn + 2t in )/a 2 is a dimensionless param- 
eter. 

We also numerically solve Eq. (|8) with V 2 = V 2 + V 2 . 
We calculate the efficiency for different values of a, tn 
etc. The results are shown in Fig. [4] (points are numerical 
results, and curve is Eq. (21)). We plot the efficiency e± 
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FIG. 4. Efficiency decay with respect to r± for transve rse d iffusion, 
the points are numerical results, and the curve is Eq. pi) . 

with respect to the rescaled dimensionless parameter r± , 
so all the points with different parameters collapse on a 
single curve. 

C. Total diffusion 

Experimentally, longitude and transverse diffusion co- 
exist during the whole process. Combining all the dif- 
fusive contributions mentioned above, we get the output 
field as (Appendix B) 

fout(k x ,k v ,t H + t) =e =^^-^-^e- Dt ^-^ 2 

(22) 

We consider input pulse with both Gaussian temporal 
and transverse profile as above, typically, Drft^ jktp are 
very small. Then the total efficiency will be 

P — \ g— 2r w g~^ T H - 

tot Y l/a H + 2/a w - 2 1 + t j_ (23) 

+ 0[(D V H 3 p , lk t p ) 2 } 

Typically, an — 1, &w — 1, so we have 

£tot — £w x £h x e R x £_L- 



D. Efficiency optimization and estimation 

Our model did not examine other decoherence pro- 
cesses, such as control field-induced scattering and 
ground state decoherence. Our results simply quantify 
the effects of motional diffusion on GEM efficiency, and 
therefore the represent upper estimates for the perfor- 
mance of GEM. 
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FIG. 5. (a) The gradient is turned off during hold time, and 
flipped for read-out. (b) An illustration of the spatial frequency 
for (7i2 (z, t), k decreases or increases depending on the sign of the 
gradient. 




t+t (S) x ID" 6 

FIG. 6. Numerical results of |<7i2(fc, t)\ with write-in gradient r\' < 
0. 

Structures with larger spatial frequency will decay 
faster under diffusion. For the ID model during hold 
time, we have (see the Appendix A) 

*i 2 (M)oc£ n (^^) (24) 
V 

The input pulse is centered at — then we have k ~ 
kjj — ki — rjti n . We see that k will increase or decrease as 
ti n increases, depending on the sign of rj. During the hold 
time [0, tn], the gradient is turned off, and k will hold its 
value. The read-out process is symmetric to the write-in 
process (see fig. [5|, returning the quasi- momentum to its 
original distribution. Fig. [6] shows the numerical result 
by solving Maxwell- B loch equations. 

One way to reduce the effects of diffusion is to remove 
the gradient during the storage part of the process, and 
only turn on the flipped gradient during readout. For 
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realistic t in and 77, we can get zero spatial frequency 
kn = 0, for which r# = 0, and this will minimise the 
diffusional decay rate. We then have Sh — an d 
Ew = e~ 2Dkitin / 3 . Including transverse diffusion, we get 
the total efficiency for input field with transverse Gaus- 
sian profile 



£tot 



-4Dk^t in /3 



1 



4D(t H + 2t in )/a 2 



(25) 



The efficiency can be improved further by choosing a 
larger transverse width a, i.e, the effects of transverse 
diffusion will be reduced by using a smooth field in the 
transverse direction. 

We note that the circumstances in which a GEM will 
be useful are those for which all dephasing, including 
that due to diffusion, is small. In this limit, a useful 
approximate expression for the GEM efficiency is given 

by 



£tot — 1 



DtHTj tp 4D(t H + 2U 



(26) 



as the inefficiencies arising from each diffusive process 
considered above add together. 

Experimental considerations give estimates of the 
achievable GEM efficiency. In particular, to ensure the 
bandwidth of the memory is large enough to absorb the 
input field, we require \rjt p \ > ^, and U n > t p to ensure 
that the whole pulse enters the medium during the write- 



in process, also 
So 



ki\ > j is required to satisfy ku — 0. 



£tot ^ 1 " 



< 1 



3 

ADt p 
3L 2 2L 2 



_ DtH _ *D(t H + 2t p ) 

2L 2 a 2 
Dt H *D(t H + 2t p ) 



(27) 



This gives a reasonable upper bound on the GEM effi- 
ciency, given the pulse duration 2t p , the hold time tn, 
and the vapour length L and beam width a. 

Experimental considerations: In experiments reported 
in p~4l [21], Rb 87 atoms were used. Typical system pa- 
rameters are ujo = 2tt • 377.10746 THz, ujo — uj c = 2tt ■ 6.8 
GHz, A = -2tt • 1.5 GHz, ^ c ~ 2tt • 20 MHz, g ~ 2tt • 4.5 
Hz, t v = 1/is, a ~ 1.45 mm, 2L = 0.2 m, rj ~ —2ir • 10 
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MHz/m, N ~ 0.5 x 10 18 m" 3 [H[2TJ[22]. The optical 
depth |/3 1 ^ 3.8 is sufficiently large. According to the 
formula in [23], we have D ~ 0.004 m 2 /s for Rubidium 
atoms in buffer gas [I4j (2TJ • 

With these parameters, the diffusive decay will be 
dominated by transverse diffusion. For example, for 
ti n = 5/iS and tn = 0, the maximum achievable efficiency 
is e tot ~ 93% (e± ~ et t, £^ = 1, ^ - 1). 

We examine the input of the (1,1) Hermite- Gaussian 
mode f in (x,y,t) ex xye~( x +y ^ a as an example of a 
higher order Hermite-Gaussian mode transverse profile. 
From Eqs. (18, 20|), we have e±_ = ( i+ T± ) 3 ? an d the longi- 



tude diffusion enects are the same as the Gaussian profile 



(i.e. the (0,0) Hermite-Gaussian mode). Thus, for diffu- 
sive decays, we have 



g (n) 
£(00) 



1 



-? 



(28) 



where is the read-out efficiency for (ij) Hermite- 

Gaussian mode. We find that the efficiency decays faster 
for higher order modes, and the ratio Eq. ([28] ) decreases 
when the storage time increases. This is in agreement 
with experimental investigations [21]. 



E. Output beam width 

After some storage time, transverse diffusion will tend 
to smear the spin wave density in the radial direction. 
Intuitively, we would expect this to lead to a spatially 
wider output beam than would be the case in the absence 
of diffusion. 

This is certainly the case when the control field is ra- 
dially uniform. To see this, we define the intensity dis- 
tribution for the read-out signal as 



I(r±) 



\fout(r±,t H + t)\ 2 dt. 



(29) 



We suppose that the control field is turned off during the 
hold time, [0, tjj] to avoid control field-induced scatter- 
ing, and that the gradient is always on and flipped at 
t = OMh- Also for typical experimental parameters, the 
effects of longitudinal diffusion is very weak, so we fo- 
cus on transverse diffusion. We solve the Maxwell-Bloch 
equation using the same method as before. For a signal 
with a Gaussian transverse profile, we find that 



I(r±) ex e 



rl/[a 2 +4D(2t in +t H )] 



(30) 



with r\ = x 2 + y 2 . Defining w r± as the width of the 



output field, we have 



D(2t in + t H ), 



(31) 



which increases linearly with storage time (Fig. [9|, at a 
rate determined by the diffusion coefficient. 

Somewhat surprisingly, the experimentally measured 
rate of expansion of the read-out signal is smaller than 
that expected from atomic diffusion by a factor of 2 to 
3 [21]. One possible explanation for this is the signal 
diffraction as suggested in [2Tj , diffusion leads to a beam 
with reduced divergence and the measurement is taken 
downstream. However in this experiment the scale of 
experimental setup is much smaller than the Rayleigh 
range, so the diffraction effect is too small to explain the 
observed discrepancy. 

Instead, we find that the anomalously narrow output 
beam width can be explained by considering the con- 
trol field with realistic transverse Gaussian profile. This 
leads to a transverse variation in the phase of the spin 
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FIG. 7. The extra phase for control field with Gaussian profile. 
Points are numerical results using typical parameters given i n th e 
main text, and the curve is the approximate expression in Eq. (|32|>. 




FIG. 8. (Color online) The intensity distribution for the read-out 
signal with tn = 16 (is. To see the expansion clearly, we have 
renormalized the maximum of I(r±) to 1, and Iq is the renormal- 
ized intensity distribution. The black solid curve is input signal, 
the blue dotted one is the read-out signal for homogeneous control 
field, and the red dashed is read-out signal for control field with 
Gaussian profile. 



wave, which, under the influence of diffusion leads to 
lower emission efficiency in the wings of the spin wave. 

To analyse this effect quantitatively, we consider 
a Gaussian transverse variation in the control field, 
f£ c (r_L) = Qe~ r - L ^ Wc , with beam waist w c . Then the two- 
photon detuning, 5, and the optical depth becomes r± 
dependent. 

From our solution for the spin wave [see Appendix A, 
Eq. (37)], we find that the inhomogeneity of the control 
field intensity will introduce a transverse variation in (5, 
which leads to a transverse dependence in the phase of 
the spin wave. Likewise, the transverse variation in the 
optical depth leads to a radially-dependent longitudinal 
shift in the spin wave 012(1", t). In combination, these 



kT" 5.8 



* 5.6 



5.4 




0.5 



2t. +t u (s) 
in H v ' 



1.5 



x 10 



FIG. 9. Expansion of the read-out signal. Circles are numerical 
results for homogeneous control field. Squares are numerical res ults 
for control field with Gaussian profile, the solid line is Eq. ( |31| for 
homogeneous control field. 



give rise to a radially dependent phase on the spin wave, 
^o(r ± ) ^ w jth the effect of the control field typically being 
dominant. We compare solutions for this inhomogeneous 
control field with solutions for the homogeneous control 
field to obtain the phase difference 6(r_\_) during hold 
time. Typically, the width of the control field, w ci is 
much larger than the width of the signal field, a, so 0(r±) 
is approximately quadratic in r±/w c : 



0(r ± ) 



21^ ti r , 



2/3 1 



2/3 In 

n 2 \ 
a v l) 



rjLU 



TjLtin + /3 



(32) 



wf. 



where f3 is the optical depth corresponding to Q. Be- 
cause of this quadratic phase variation across the spin 
wave, diffusion acts to wash out the spin- wave coherence 
more quickly at larger radius, so the read-out efficiency 
is suppressed at larger r±. This will tend to reduce the 
apparent width of the emitted read-out signal. 

Experimental considerations: In the experimental re- 
sults reported in [21], w c ~ 3 mm, and U n ~ 2 /xs. Using 
these parameters, Fig. [7] shows the transverse variation 
in the phase of the spin-wave at {z = 0, t = OMh) in 
the absence of diffusion. When diffusion is introduced, 
this transverse phase variation is smeared out, leading 
to reduced read-out efficiency in the wings of the spin- 
wave. Figure [8] compares the numerical results for the ex- 
pansion of the read-out signal after a specific hold time, 
tu = 16 /is, with a homogeneous control field (dotted, 
blue) and with a spatially varying control field (dashed, 
red), assuming the diffusion rate D = 0.004 m 2 /s. Fig- 
ure [9] shows the variation in the width of the output field 
as a function of hold time. We see that the expansion is 
slowed for a control field with Gaussian profile (squares) , 
compared to the case of a uniform control field (circles). 
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Importantly, this corresponds to a reduction of the beam 
width expansion-rate by a factor of 2. The apparent dif- 
fusion rate extracted from this slower expansion rate is 
D e ff ~ 0.002 m 2 /s. This is quantitatively in agreement 
with the observations in [21] . 

V. SUMMARY 

We have studied the effects of diffusion on the efficiency 
of the A-gradient echo memory, both numerically and an- 
alytically. We find that the efficiency is dependent on 
the spatial frequencies k for both longitude diffusion and 
transverse diffusion: higher k leads to more pronounced 
diffusive effects, and reduced efficiency, as expected. We 
show that the storage efficiency can be improved by ap- 
propriate choice of the gradient during the hold phase. 

We established a mechanism by which the rate of ex- 
pansion of the transverse width of the beam is reduced, 
compared to the naive expectation of diffusive effects. 
This mechanism arises from the effects of diffusion on the 
transverse variation in the spin wave phase. We showed 
that with an experimentally reasonable choice of param- 



eters, the magnitude of this effect is the same as that 
observed in recent experiments. When the density of the 
buffer gas in increased, the collision rate increases, lead- 
ing to a smaller diffusion rate. However, this will lead to 
collision-induced dephasing, which will dominate at suffi- 
ciently high buffer gas pressures. This implies a trade off 
between diffusion- and collision-induced dephasing. This 
will be the subject of future research. 
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APPENDIX A 

The Maxwell-Bloch equation for the 1-dimensional 
model is 



E(z,t) 

i5(z,t)a 12 (z,t) + DV 2 z a 12 (z,t), 



+ i^-E(z,t). 



(33) 



cA 

To find the solution during [— to,0], we first solve the 
equation without diffusion, then introduce the diffusion 
effects to our solutions. 

When D = 0, we can make transformation 



a 12 (z,t) =e % ^ "e 



< k °- k ^ z a 12 (z,t), 
E(z,t) =e- l ^ z E(z,t), 
and get the new equations 



(34) 



3 z E{z,t) =i^£l^a 12 (z,t), 
d t a 12 (z,t) = - irjza 12 (z,t) + ig eff E(z,t) 



(35) 



where g e ff — gQ c / A. Following the method given in [TO] , 
and using the boundary conditions a 12 (2;,£ — >• — oo) = 
and E(z = —L,t < 0) = fi n (i), we integrate the first 
equation and substitute it in the second one. Making use 
of Fourier transformation, we find 

E(k,t) = J in {- + A r +t)\-\-^- 1 G( V ,/3,L), (36) 

T] T]L T] 



and 



G%frL) = ^e-l^/ 2 S i2iii(7r^|)|77L|-^r(i/3), 



where E(k,t) = J E(z,t)e~ ikz dz, f3 



jeff 



N 



is the op- 



tical depth and we assume /3 is sufficiently large, T(if3) 

is the Gamma Function, fi n (t) = fi n (t) elE ^ L is the in- 
put pulse. According to the Maxwell-Bloch equations, 
we have a 12 (k,i) = -*j^E{k,t). 



eff 



We transform <r 12 (A;,£) back to cr 12 (/c,t), 
, k kj 



a 12 (k : t) =f in (- 



V 



V 



7]L T] 



vL'g ef fN 

(37) 



with ki = i£ 



'0 — K c ~ l • 

Now we introduce the diffusion, for the short time in- 

Dk 2 At 



k - k, 
uce th« 

terval [t,t + At], diffusion will cause a decay e 



to 
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a(k,t), or equally, there will be a decay e ~ Dk2At on the 
signal fin(t') with k = ki — rj(t — t'). So the total decay 
during the write-in process for fi n (t') is 



-D f0(k i -r l (t-t , )) 2 dt 



e 3r ? 



^(fcf-^+ryt') 3 ) 



Thus, the solution for ayi at t = is 



/, _ \ ~ D (h 3 h 3 \ „ ,K — Ki. a9 Z N 

en (M) =e^r (fc * fe >/ in ( e 1 ^ 

77 



?7 7?L 



zf3 sgn( 



k ki (3 ^ c q 



7] r]L g eff N 



(38) 



We have assumed that the bandwidth of the memory 
is larger than the bandwidth of the input signal, \rjL\ ^> 
Au s , and the optical depth is sufficient large, f3 > 1. 
The signal will be absorbed near z = 0, and 012(2,0) 
and E(z,0) is nonzero only near z = 0, so we can treat L 
as infinity during [0, Also notice that, the gradient 
is turned off during [0,£#], and the spatial frequency k 
will hold its value. 



To get the solution in [0, we solve Eq. (33) with 
initial condition <Ji 2 (k,t = 0) for 012 and open boundary 
condition for E. In k space, we find 



012 (Mif) = e 



-Dk z tJ 



'elf 1 



tH <T 12 (k,0) 



(39) 



where k = + ko — k c . Notice that 012 (&,£#) get a 



phase e 
v g (k) = 



cA 

C h-h. H 



< ff N r 

c(k—k) : 



, so the group velocity for <Ji 2 (z,t) is 
. If the memory broadening \rjL\ is not 



much larger than the signal pulse bandwidth, the spin 
wave cr 12 (z,t) will be nonzero near the ensemble bound- 
ary. Then the spin wave will propagate to the boundary 
and be reflected, this may ruin the spin wave coherence 
near the boundary and lower the memory efficiency. One 
way to avoid this effect is turning off the control field dur- 
ing storage, which makes the effective coupling g e ff = 0, 
and the group velocity v g = 0. 

To find the values for 012 and E in the duration 
\Ph^h + to], one needs to solve a modified version of 
Eq. ( [33] ) where the sign of irjz is reversed. We fol- 
low the method given in [10 , propagate these equation 
backwards with final conditions E(z = L,t > tn) = 
i out (t) •> 12 (z,t 00) = 0. Similar to the write-in pro- 
cess, at time we have 



4 3 N T 



D (h 3 h 3 ^ n / k — ko. 

-V 

x k-h _ A|-i/5 sgn( *^*L _ A)_^ G * 

T] T]L 7] T]L g e ffN 

(40) 



By matching the two solutions for 012 at tn Eqs. (|39j), 
(Eol), we get 



where 
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a = |/,/. ( / - ^ » i-' 2 ^> i 



»ef f N 

e- c < t+ £) Y{iP)/Y{-ifi 



is a phase factor, d w (t) = e 3 ° ^ vtj3 \ d H = 

z-Dik^tftn and dR ^ = ^(kz-ikt-mf) are the diffu _ 

sion decays for the write-in process [—to, 0], storage time 
[0, and read-out process \pH^H + ^0] respectively. 



APPENDIX B 

The Maxwell-Bloch equation for the 3-dimensional 
model is 

a 12 (r,t) =i 9 -^e^- k ^E{v,t) 

- (ir ? z)<j 12 (r,t) + J DV 2 (Ji2(r,t), 

To solve these equations, we first transform transverse 
coordinates x, y to Fourier space k x ,k y , 

(r 12 {k x ,k y ,z,t) = - (irjz + j k )a 12 (k x ,k y ,z,t) 

H~ i—^—e'^ c ^ E{k x ,ky,z,t} 

+ D\7 2 z a 12 (k x ,k y ,z,t), 



dz 



E(k x ,k y ,z,t) =i — r-^e 



9 NSl c r -i(kp-k c )z 

cA 



012(kx,ky,Z,t) 



.g 2 N 



+ i E(k x , k y ,z,t). 



(43) 

where 7^ = D(k 2 + ky). Now we make the following 
transformation: 

0i2fe, k y , z, t) =e lkt (T 12 (k x , k y , z, t), 
E(k x , ky, z, t) =e^ ki E(k x , ky, z, t), 
then we have 

(?\2(k x , k y , z,t) = - (if]z)a 12 (k x , k y , z, t) 
.gSl 



(44) 



+ i^e i fo- k J*E(k x ,k v ,z,t) 
+ DV 2 z a 12 (k x ,k y ,z,t), 



l z E(k x ,k y ,z,t) =i 9j^L e -i( k0 ->. c y, 



(Ji 2 (k x ,k y ,z,t) 



,g 2 N - 
+ i E(k x , k y , z,t). 



(45) 



These are actually quasi-lD equations, so we can solve 
these equations by the method we used before, and the 
output field is: 



fout(t H + t) =d w (t)d H d R (t)f in (-t)G (41) f out (k x ,k y ,t H + t) = d w (t)d H d R (t)f in (k x ,k y ,-t)G. 
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We transform back to f ou t(k x , k y ,tu + t), and get 
fout(k x ,k y ,t H +i) =d w {t)d H d R {t) 

(46) 

xd ± (t)f in (k x ,k y ,-t)G, 

where d±(i) = e~ 2lkt e~ lktH is the transverse difusion 
decay. 



